THE CRYSTALLOGRAPHY OF MARTENSITE TRANSFORMATIONS 
III. FACE-CENTRED CUBIC TO BODY-CENTRED TETRAGONAL 

TRANSFORMATIONS* 

J. S. BOWLES and J. K. MACKENZIEf 



The theory developed in Parts I and II is applied to transformations from face-centred cubic to 
body-centred tetragonal in the alloys Fe-C, Fe-Ni, Fe-Ni-C. The theoretical predictions of habit 
planes, orientation relationships and directions of the homogeneous strain are found to be in satisfac- 
tory agreement with all the available experimental data. Thus, the hypothesis that the inhomogeneous 
shear is a part of the twinning shear in the final structure has been substantiated for the transforma- 
tions considered. 

The habit plane interface envisaged by Frank is shown to be a special case of the type of interface 
expected on the present theory : all planes in the zone of the invariant line, and in particular the 
twinning plane, meet edge to edge in the interface. 

LA CRISTALLOGRAPHIE DES TRANSFORMATIONS MARTENSITIQUES III. LES 
TRANSFORMATIONS DU SYST£ME CUBIQUE A FACES CENTR£ES EN TETRAGONAL 

CENTRE 

La theorie developpee dans les Parties I et II est appliquee aux transformations du systfeme cubique 
& faces centrees en tetragonal centre, dans les alliages Fe-C, Fe-Ni, Fe-Ni-C. II est constate que les 
predictions theoriques des plans Iimites, des relations d’orientation et des directions de deformation 
hoinoyhnc sont en accord satisfaisant avec toutes les donn6es experimentales, disponibles. Ainsi, 
l’hypothfese que la tension non-homog£ne fait partie de la tension de maclage dans la structure finale, 
a 6te confirmee pour le cas des transformations considerees. II est montre que l'interface constitute 
par un plan limite, envisagee par Frank, n’est qu'un cas particulier de l’interface proposee par la 
presente theorie, k savoir: tous les plans dans la zone de la ligne invariante, et en particulier le plan 
de maclage, se rencontrent suivant une arete commune dans l’interface. 

DIE KRISTALLOGRAPHIE DER MARTENSIT TRANSFORMATIONEN. III. TRANS- 
FORMATIONEN VON KUBISCH-FLACHENZENTRIERTEN ZU TETRAGONAL- 
RAUMZENTRIERTEN GITTERN 

Die in Teil I und II entwickelte Theorie wird auf die Umwandlungen kubisch-flachenzentriert 
zu tetragonal-raumzentriert in Fe-C, Fe-Ni und Fe-Ni-C Legierungen angewandt. Die theoretischen 
Voraussagen fiber Habitusebenen, Beziehungen der kristallographischen Orientierungen und Rich- 
tungen homogener Verzerrungen stimmten zufriedenstellend mit alien vorhandenen experimentellen 
Angaben Uberein. Somit ist die Hypothese, dass die inhomogene Scherung einen Teil der Zwillings- 
scherung in der Endstruktur darstellt, fur die hier betrachteten Transforma tionen bestarkt worden. 

Es wird gezeigt, dass die von Franck vorgeschlagene Habitusebenengrenzflache ein Spezialfall des 
nach der vorstehenden Theorie zu erwartenden Grenzflachentyps ist. Alle Ebenen in der Zone der 
invarianten Richtung und im besonderen die Zwillingsebene schneiden sich in der Grenzflache. 



1. Introduction 

In this paper the geometrical theory of marten- 
site transformations which was developed in 
Parts I and IIJ [1] will be applied to transforma- 
tions from face-centred cubic to body-centred 
tetragonal lattices, and its predictions of habit 
planes, orientation relationships and other geo- 
metrical features will be compared with the 
available experimental data for iron -carbon, iron- 
nickel and iron-nickel-carbon alloys. 

In I and II a phenomenological theory of the 
relations between the various crystallographic 
features of martensite transformations was devel- 
oped. This theory aims at describing the total 
atomic displacements produced by martensite 
transformations but not necessarily the motions of 
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JThese papers will be referred to as I and II. 



atoms during transformation. It has been shown 
that the total atomic displacements can be 
described consistently by means of a homogeneous 
strain which is composed of an invariant-plane 
strain on the habit plane and a small dilatation, 
followed by an inhomogeneous shear which pro- 
duces no further observable change in shape. 
This shear can only be homogeneous within small 
volumes and it was proposed that within such 
volumes it is a part of a twinning shear in the final 
lattice. It has been shown that these conditions, 
together with a knowledge of the correspondence 
between the two lattices, their dimensions, and the 
twinning plane involved, is sufficient to determine 
the total strain in terms of a single parameter, 
6 , and, further, that the resolution of this total 
strain into component strains of the required 
types is uniquely determined. 

In §§2, 3, of the present paper the correspon- 
dence between the initial and final lattices, and 
also the relevant twinning plane, are shown to be 
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the same for all the transformations considered. 
In these sections the experimental origins of this 
information are described and quantities which 
are required for direct substitution into the 
equations of the theory are derived. In §4 the 
parameter 6 is introduced and substitution into 
the main equations of the general theory is carried 
out. To apply the theory to a particular transfor- 
mation it is necessary to find the appropriate 
value of 0 by calculating the variation of the habit 
plane with 6 and determining that value for which 
the predicted and observed habit planes agree. 
Of course the extent to which there is agreement 
between the observed and the theoretically pos- 
sible habit planes is a first test of the hypothesis 
that the inhomogeneous shear is a part of a twin- 
ning shear. The actual comparisons of predicted 
and observed habit planes are carried out in §5 
and the appropriate values of 6 are determined for 
the various transformations. The remaining sec- 
tions of the paper are concerned with detailed 
analyses of individual transformations and show 
that the predictions of the theory are in substantial 
agreement with experiment. 

The notation used in this paper is the same as 
that used in I and II. 



2. The Correspondence and the Principal Axes 



Jaswon and Wheeler [2] have determined the 
correspondence between f.c.c. and b.c.t. lattices 
for those transformations in iron-carbon alloys 
which produce the Kurdjumow-Sachs orientation 
relationship [3]. For the variant of this relationship 
shown in Figure 1, their correspondence is 



(2.1) 



(bCf) = 



0 

0 

1 




where the bases b and f define the b.c.t. and f.c.c. 
unit cells respectively; the corresponding metrics* 
are (b*Gb) = (diag: a' 2 , a' 2 , y 2 a' 2 ) and (f*Gf) = 
(diag: a 2 , a 2 , a 2 ) where a' and a are the respective 
lattice parameters and y is the axial ratio of the 
tetragonal cell. This correspondence was chosen 
by Jaswon and Wheeler because it ensures that 
the atomic displacements involved in the trans- 
formation are as small as possible. It has been 
shown [4; 2] that the arrangement of interstitial 
carbon atoms in iron-carbon martensite is precisely 
that expected from this correspondence. It has 
also been shown [5] that the same correspondence 

*The notation (diag: a, b, c) represents a square diagonal 
matrix with elements a, b, c in the leading diagonal, all other 
elements being zero. 



results from a consideration of the analogous 
variants of all orientation relationships between 
Nishiyama [6] and Kurdjumow-Sachs, including 
that of Greninger and Troiano [7]. 




Figure 1. Stereographic projection showing the variant 
of the Kurdjumow-Sachs orientation relationship considered 
in this paper. The relevant twinning plane and twinning 
direction are also shown. The subscripts F and B refer to the 
f.c.c. and b.c.t. lattices respectively. 



The correspondence and the metrics of the two 
bases b and f determine the principal axes, and 
the associated principal strains of the total strain, 
as described in 1 1- § 4. In the present case the 
results are easily obtained by inspection, for the 
mutually orthogonal directions [101] F , [10l] F , 
[010] p correspond to the mutually orthogonal 
directions [100] B , [010] B , [001] B . The former set of 
directions can therefore be chosen as principal 
axes and the basis p, comprised of unit vectors in 
these directions (Figure 2) defined to be [f; pi] = 




Figure 2. Stereographic projection showing the orienta- 
tion of the bases P, I, L, K relative to the f.c.c. lattice. 
The invariant line and the invariant normal make angles a 
and (3 respectively with [001 ] F . 




226 



ACTA METALLURGICA, VOL. 2, 1954 



[1/V2, 0, 1/V2]fi [f; p 2 ) = [1/V2, 0, -1/V2 ]f, 

If; pa] = [0, 1, 0] K . 

Thus, 

/1/V2, 0, 1/V2\ 

(2.2) (pTf) = 1/V2, 0, -1/V2 = (fTp)\ 

V 0 , 1, 0 / 

The elements of the diagonal matrix M, represent- 
ing the strain referred to the basis p which extends 
the base vectors p to their final lengths without 
rotation, are simply the ratios of the final to the 
initial identity periods along the principal axes. 
Thus, 

(2.3) M 2 = {a' /a) 1 (diag: 2, 2, y 1 ). 

3. The Twinning Plane and Direction 

Greninger and Troiano [7] have observed 
twinned martensite plates in an iron-22% nickel- 
0.8% carbon alloy and have determined the 
indices of the twinning plane. Examination of 
their pole figure shows that the operative twinning 
plane is (112) B for the variant being considered 
in this paper. 

That the twinning plane is the same for trans- 
formations which lead to the Kurdjumow-Sachs 
orientation relationship can be shown as follows. 
This orientation has the property that its variants 
are twin-related in pairs. The twinning shear 
relating such a pair variants can be obtained by 
stereographic analysis [5] making use of the 
correspondences for these variants. The relevant 
twinning plane is again found to be (112) B and the 
associated twinning direction is [llI] B - This twin- 
ning direction is the same as that obtained from 
the general equation for the twinning shear, 
II-(7.7), which gives 

(3.1) (bS t b) = I - f[\, 1, 1](l, 1, 2), 

where / = (2 — y i )/2{y- + 2). 

The intermediate basis I used in the general 
theory can now be determined explicitly. The 
base vector i 2 ( = u) is a unit vector parallel to the 
direction in the initial lattice from which the 
twinning direction is generated. The base vector 
i 3 (= h) is a unit vector normal to the plane from 
which the twinning plane is generated, while ii 
is a unit vector perpendicular to both. Using the 
correspondence (2.1) it follows that the twinning 
direction [llI] B is generated from the direction 
(bCf^HHIIb = [2 2 0] F , so that [f ; i 2 ] = [l/v% 
— 1/V2, 0] F . Similarly, the twinning plane (112) B 
is generated from the plane (112) b (bCf) — (220) F 



so that [f; i 3 ] = [1/V2, 1/V2, 0] F and [f; ij] = 
[0, 0, 1] F . These base vectors are shown in Figure 2. 
Finally, using (2.2) to transform the components 
of these base vectors to the basis p it follows that 
the transformation 

/ 1/V2, 1/2, 1/2 \ 

(3.2) (pTi) = I — 1/V2, 1/2, 1/2 1 = (iTp)'. 

\ 0 , — 1/V2, 1/V2/ 

It should be noted that the twinning plane 
(112 ) b is generated from a plane of symmetry 
(110) F as required by the general theory, II- § 7. 

4. The Invariant Line Strain 

The invariant line strain S that is obtained from 
the total strain S< by removing a suitable dilatation 
1/5, is calculated in II-§§6, 9 and is expressed 
ultimately in terms of the elements g r „ of a matrix. 

, (4.1) G = (iTp) 5 2 M 2 (pTi). 

The dilatation 5 is the only unknown quantity in 
this equation, and is therefore the only adjustable 
parameter in the theory. However, 5 itself is not 
the most convenient parameter to use. Equation 

(2.3) shows that 5 2 M 2 = 6 2 (a'/a) 2 Mo 2 , where Mo 2 = 
(diag: 2, 2, -y 2 ), which for a fixed y is independent 
of the sizes of the initial and the final cells. Thus, 
if a new parameter 0 is defined by the relation 

(4.2) 0 = 5 (a'/a), 

all transformations with the same y can be treated 
simultaneously for they differ only with respect 
to the value of 0. In terms of the parameter 0 

/ 2, 0 ,0 \ 

(4.3) G = (g„) = 0 2 0, 1 + 7 2 /2, 1 - 7 2 /2 , 

Vo, 1 - 7 2 /2, 1 + 772/ 

(4.4) 

* i A- 0 - 0 \ 

G 1 = (gr.) = o, (2 + 7 2 )/47 2 , (t 2 — 2)/47 2 ), 
\0, (7 2 - 2)/4 7 2 , (2 + 7 2 )/4 7 V 

and 

(4.5) |G| = g = 47 V. 

In II it was shown that the matrix representing 
the invariant line strain S has a simple form when 
referred to a basis L, such that h is parallel to the 
invariant line. The base vectors 1 1( 1 2 , 1 3 are shown 
in Figure 2. The transformation 

( cos a, — sin a, 0\ 
sin a, cos a, 0 I, 

0 , 0,1/ 
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where a(|a| < 90°) is the angle between li and 
[001] F = ii, and is determined by equation II- 
(6.10). On substituting the g rt given in (4.3) into 
this equation it is found that 

(4.7) cos a = (1 - 0 2 - e\ 2 /2) h /e{\ - y/2)\ 

(4.8) sin a = ( 26 2 - 1)*/0(1 - y72)*, 

giving two values of a which are merely opposite 
in sign. It should be noted that these equations 
only determine real values of a provided 1 / \/2 < 
0 < 1/(1 + 7 2 / 2 )*, thus imposing limits on the 
possible values of 0; the minimum and maximum 
values of 0 correspond to a = 0° and 90° respec- 
tively. 

The matrix representing S referred to the basis 
L is given by II-(8.7), (6.18) 

(4.9) (lSl) = 

( 1, (g22 — gu) COS a sin a, 

0, gg *3 COS CD 

0, gg*f sin 03 

g 23 Sin a \ 

— (gg*l cos a cos « + sin «)/gtl 1 , 
(- gg*l cos a sin w + cos w)/gtl/ 

where o> is the angle between the initial and final 
positions of the twinning plane and is found 
from II-(8.5) to be 

(4.10) cos 03 — (207 + 1 + 7 2 /2)/(l + 07 )( 2 + 7 1 )*, 

giving two values of <o, opposite in sign. It should 
be noted that cos 03 must be positive since the 
criterion for selecting the correspondence ensures 
that the total strain is as small as possible: if cos w 
were negative the twinning plane would have to be 
rotated by more than 90° by the transformation. 
Considerations such as this show that all the 
square roots must be taken positively unless the 
contrary is stated. 

The two values of a taken together with the two 
for co determine four possible invariant line strains, 
all of which have the same principal axes and 
associated principal strains, and the twinning 
plane is generated from the same plane (110) F 
in all four cases. The solution obtained by taking 
a and co both positive, say, will be referred to as 
the (a+, co +) -solution and so on. In II-§ 8, it was 
shown that the (« + , co — )-solution refers to the 
twin of that variant of the (a + , co+) -solution 
which is obtained by a rotation of 180° about the 
normal to the (110) F symmetry plane. The (a-, 
co-)- and the (a — , co+)-solutions are related in 



the same way. In general, the (a + , co-j-)- and 
(a — , co +) -solutions are distinct, but in the present 
case they are variants related by a rotation 
through 180° about the normal to the (001) F 
symmetry plane. The other pair are, of course, 
related in the same way. Examination of Figure 2 
shows that this symmetry operation leaves the 
principal axes, principal strains and the plane 
(110) F invariant* while it changes 1 1 (a + ) into 
li(a — ). This result can be demonstrated formally 
by showing that the (a—, co+)-solution, when 
referred to the variant of its L basis, is identical 
with the (a+, co+)-solution. 

The factorization of S can now be carried out 
making use of the relations given in II-§ 9. In 
the present class of transformations certain simpli- 
fications arise and the formulae can be given 
geometrical interpretations which are sometimes 
useful. For the (a+, co+)-solution, unit vectors 
along the normal to the habit plane and in the 
direction of the displacement in the invariant 
plane strain are given by 

(Pi L) = (0, -1, cos /3) l /(1 + cos 2 (8)*, | 

C 4 - 11 ) (p ; ^ _ (sin^ — Qog^cog^j /(l-f-cos 2 /3)*,j 

[k; d] = [0, cos a, -1] K /(1 + cos 2 a) 1 , 

(4.12) jj. _ [gin l 8 j cosci!, — cos /3]i /(l+cos 2 tt)*, 

where k is the basis introduced at the end of 
II-§6 and 0, the angle between [100] 1 and the 
invariant normal of S, is given by II-(6.11) 

(4.13) sin 8 = V2 t( 20 2 - l) J /(2 - 7 2 )*- 

It follows from (4.11) that for constant 8 the habit 
plane is a fixed plane in the L-basis. Since the 
L-basis only varies by a rotation about [110 ]f, 
the normal to the habit plane varies in the same 
way for constant Similarly for constant a the 
direction only varies by a rotation about [110] F . 
Further, the magnitude m (I-§5) of the invariant 
plane strain reduces to 

(4.14) m = V 2 0(1 - 07), 

while the fraction k of the twinning shear involved 
in the complementary strain is given by 

(4.15) 1 + 2 k = cos a cos p. 



‘Since [001] F is a line of symmetry in the initial position of 
the twinning plane (110) F , these facts show that the two 
possible values of a can differ only in sign. 
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5. The Habit Plane 

It was shown in II-(9.3) that the habit plane 
referred to the basis l is of the form (0, sin u, 
gu*"^ — cos oj) l . On substituting for w and g 33 * 
this becomes 

(5.1) (0, sgn (1 - 772)*, - [I-tW-D/^*)!,. 

The habit plane can now be referred to the basis F 
using the transformation 

(fTl) = (fTp)(pTi)(iTl), 

(5.2) / sina/V2, cos a/y/2, 1/V2\ 

= 1 — sin a/\/2, — cosa/VS, 1/V2 1 . 

\ cos a —sin a, 0 / 

The variation of the habit plane with 9 for 7 = 1, 
i.e., for body-centred cubic structures, is shown 
in Figure 3 for all four solutions, and the relation- 



a=0* 




Figure 3. Stereographic projection showing the variation 
of the habit plane with 6 for 7 = 1. As 9 1 varies between the 
limits 1/2 and 2/3, a varies from 0 to 90 degrees. The open 
circles represent the four solutions obtained for the habit 
plane for some particular value of 9. 

ship between the solutions is clearly that described 
previously. Since the four habit planes are crys- 
tallographically equivalent, it is only necessary to 
consider one eighth of the whole projection in 
order to compare experimental habit planes with 
those predicted by the theory. Accordingly, in 
Figure 4 the variation of the habit plane with 6 
and 7 has been plotted for the (a+, w+) -solution. 

The compositions of those alloys for which 
habit planes have been determined are given in 
Table I, which also shows the corresponding values 
of 7 and ( a' /a ) 2 . In all the experimental work, 
with one exception [7] only the general indices 
of the habit plane have been determined. In 




Figure 4. Stereographic projection showing the variation 
of the habit plane with 9 and 7 . 

Figures 5, 6, 7, where the measured habit planes 
are compared with the theory, the experimental 
habit planes have been plotted in the stereo- 
graphic triangle containing the predicted curves. 
It will be shown in the next section that Greninger 
and Troiano’s habit plane does indeed lie in this 
triangle. 

It is evident from these diagrams that the 
experimental results are in general agreement with 
the theory. The observed habit planes lie close to 
the calculated curves for the relevant 7, and the 
predicted lateral shift of the habit plane curves as 
7 increases is in good agreement with observation 
(Figure 5). 




Figure 5. Comparison of the predicted habit planes with 
those measured by Greninger and Troiano [7] for various 
iron-carbon alloys. Following Greninger and Troiano, the 
results for 0.95%C. have been plotted in the adjoining 
triangle. According to the theory each group of points should 
lie somewhere on the curve for the appropriate 7 . 

The exact location of a habit plane on a particu- 
lar curve is determined by the value of the para- 
meter 9 = 5a' /a. Although the exact values of 5 
are not known, they must be close to unity (I-§4). 
Thus, by regarding 9 as equal to a’/a, a crude 
comparison of the predicted with the observed 
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effects of varying 6 can be made. The theory then 
predicts that a' /a decreases either as 7 increases 
or as the habit plane is displaced along a curve 
towards (010) F . Examination of the data in 
Table I shows that a' /a always decreases as 7 
increases. Furthermore, in those transformations 
for which 7 = 1 (Fig. 7) and also in those for which 
7 = 1.041, 1.045 (Fig. 6 ) the habit plane moves in 




Figure 6. Comparison of the predicted and measured 
habit planes for transformations with 7 ~ 1.045. 




Figure 7. Comparison of the predicted and measured 
habit planes for transformations with 7 = 1 . The results for 
30.0% and 32.5% Ni are those obtained respectively by 
Machlin and Cohen [14] and Greninger and Troiano [7], 

the predicted direction as a' /a varies. Thus, the 
theory accounts for the general trends of habit 
plane variation. Quantitative agreement is achieved 
by evaluating the adjustable parameter S. 

The values of 0 2 which lead to the best agreement 
with the experimental habit planes have been 
estimated for each transformation, and are given 
in Table I, together with the values of S derived 
therefrom. The values of 5 are all within 1.5% of 
unity, in agreement with the observational evi- 
dence that any dilatation of the habit plane could 
not exceed a few per cent. (I-§4). 

The feature of the experimental data that seems 



TABLE I 



Composition 


y 


(s'/o)’ 


0* 


b—da/a ' 


References 


Carbon Steels 
0 %C. 


1.000 


0.650 


0.6667t 


1.013 




0.35%C. 


1.016 


0.642 


0 . 6596f 


1.014 




O.02%C. 


1.041 


0.630 


0.6481f 


1.014 


Roberts* [8] 


1.40%C. 


1.063 


0.620 


0.63861 


1.015 




1.78%C. 


1.080 


0.612 


0.620 


1.006 




Fe-30% Ni 


1.000 


0.642 


0.647J 


1.004$ 


Owen, Yates and 


Fe-32.5%Ni 


1.000 


0.640 


0.645 


1.004 


Sully [9]; Bradley, 
J ayandTay!or[10] ; 


Fe-0.8%C. 


1.045 


0.627 


0.630 


1.002 


Jette and Foote 
[111; Phragmen; 
[12]; Nishiyama 
[61. 


-22% Ni 










Greninger and 

Troiano [7J. 



*For Fe-C alloys Roberts summarises all the existing data 
and gives a = (3.548 + 0.044x)kX, a ' = (2.861 - 0.013*) 
kX, 7 = 1.000 + 0.045*, where * is the wt. % carbon. For 
pure iron, Barrett [13] gives a = 3.564A = 3.557kX differing 
from Roberts’ extrapolated value and obtained by extrapola- 
tion of high-temperature data. 

t® 2 = i/a + w 2) = ew 

jThe wide scatter of the observed habit planes makes it 
difficult to estimate 0 2 . The value given has been obtained 
by assuming S = 1.004 (as for Fe-32.5% Ni) and using the 
value of (ar/a)* to calculate 8 1 . The theoretical habit plane 
for 0 2 = 0.647 lies within the scatter of the observations. 

most difficult to reconcile with the theory is the 
wide scatter of the observed habit planes, par- 
ticularly in the iron-nickel alloys (Fig. 7). Since it 
seems unlikely that the whole of this scatter can 
be due to experimental error other possibilities 
must be considered. The scatter parallel to the 
predicted curves could, of course, be attributed to 
variations in 9 from plate to plate. Such variations 
could arise from changes in 5 and/or a' /a due to 
inhomogeneneity of composition. The scatter 
normal to the predicted curves cannot be explained 
in this way, and, if it is a true effect, it may be 
necessary to abandon the hypothesis that the 
complementary shear occurs on the twinning 
plane and in the twinning direction and examine 
the consequences of fulfilling only one of these 
conditions (I-§§6, 9). However, it is clear that all 
the general features of habit plane variation are 
accounted for satisfactorily by the present hypo- 
thesis, and in the remaining sections it will be 
shown that by using the values of 9 determined 
above, other crystallographic features of the 
transformations can also be accounted for quanti- 
tatively. 

6. Greninger and Troiano’s 
(259) -Transformation 

The most complete experimental analysis of a 
martensite transformation is that made by Grenin- 
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ger and Troiano [7] of the transformation in 
Fe-0.8% C-22% Ni. They determined the habit 
plane and orientation relationship for a particular 
variant and also the direction of the homogeneous 
strain, assuming this strain to be a simple shear 
on the habit plane. The predictions of the theory 
for this transformation will now be compared with 
their observations. 

It is clear from Figure 6 that the predicted and 
experimental habit planes are in good agreement 
for 7 = 1.045, if 6 2 = 0.63. Of the four possible 
strain matrices which can now be obtained using 
these values it is only necessary to consider any 
two which characterise a parent and a variant 
of its twin; the other two solutions are just real 
variants of these. A suitable pair is the («+, o>+)- 
and the (a+, co — )-solutions, and of these it is 
the (a + , co +) -solution which leads to agreement 
with Greninger and Troiano’s orientation relation- 
ship. For this transformation a = 72° 26.87' and 
co = 4° 57.30' and the (a+, co+) -solution* is 

fl, - 0.082238, 0.272696\ 

(6.1) (lSl) = ( 0, 1.103662, 0.016254 , 

\0, 0.095683, 0.948345/ 

while the transformation (fTl) is obtained from 
equation (5.2) with sin a = 0.953443 and cos a = 
0.301574. 

The relative orientation of the two lattices can 
now be found very simply, for the indices, relative 
to the basis F, of any direction [u v w] B are (fSf) 
(fCb)[w v w] B while a plane (hkl) B has indices 
(M/) b (bCf)(fSf ) _1 relative to f. In Table II the 



TABLE II 

Comparison of Predicted and Measured Orientation 
Relationships in Fe-22% Ni-0.8% C. 



Corresponding 

directions 


Measured 

angle 


Predicted 

angle 


[1T0]f — [ llI] B 


2.5° 


2.7° 


[011]f — [U1]jb 


6.5° 


6.6° 


[121]f — [1011b 


2.0° 


1.9° 


(111)f — (101) B 


about 1° 


0.2° 



predicted angles between various corresponding 
directions and planes are compared with those 
measured by Greninger and Troiano [7 ; Figs. 3 and 
6]. It is evident that the two sets of angles are in 
very good agreement. The greatest discrepancy 

*The six-figure accuracy is needed to obtain self-consistent 
numerical results. 



(0.8°) is in the angle between (111) F and (101) B 
but this is of doubtful significance since it is not 
indicated in Greninger and Troiano’s paper whether 
the angle was measured directly or obtained by 
graphical construction. Since a gnomonic projec- 
tion was used the latter possibility seems the more 
likely, and in this case the discrepancy is not 
surprising. Notwithstanding its magnitude, the 
predicted displacement of the (101) B -pole from the 

(111) p-pole is, as far as can be judged, in the 
observed direction. The other three angles appear 
to have been rounded off and the values given are, 
in fact, inconsistent with a tetragonal structure 
with 7 = 1.045. (The results imply an angle of 
not more than 69.0° between [111] B and [llI] B 
whereas the correct angle is 69.3°). Thus, the small 
discrepancies in these three angles may be even 
further reduced, and since these three angles 
define the fourth, the estimate of 1° may need 
some amendment. 

For this transformation the theory not only 
predicts the observed orientation relationship but 
also the observed orientation of the habit plane 
relative to both lattices. It is readily demonstrated 
that a pole figure showing the austenite axes, the 
predicted martensite axes and habit plane, can be 
rotated into coincidence with the pole figure 
determined by Greninger and Troiano [7; Fig. 6]. 
The theory, therefore, predicts the correct habit 
plane: not merely a plane of the correct form. 

The orientation relationship derived from the 
(a + , co — ) -solution is simply a variant of the 

(112) B -twin of the orientation just considered. 
Greninger and Troiano have observed twins on 
this plane so that all orientations predicted by the 
theory actually occur. However, the reported 
orientation is produced in preference to its twin. 
A rational explanation of this observation becomes 
apparent on carrying out the factorization into 
component strains of both the (a+, «+)- and the 
(a+, w — ) -solutions. The homogeneous strains are, 
of course, variants of the same strain, while the 
inhomogeneous strains differ with respect to the 
fraction of the twinning shear involved. It is found 
[I I- (9.8)] that 0.4077 of the twinning shear is 
involved in the (a+, co+)-solution, while the other 
solution involves the fraction 0.5923. It is reason- 
able to suppose that the twin orientation is less 
favoured because of the larger fraction of the 
twinning shear that is involved in its formation. 

The only experimental data which remain to be 
considered concern the magnitude and direction 
of the homogeneous strain. According to the theory 
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(I I-§9) this strain is represented when referred to 
the basis F by the matrix 

(6.2) (FpiF) = I + m dp', 
where 

(6.3) -d = [0.21017, 0.61850, - 0.75715] P , 

(6.4) p' = (0.18726, 0.55102, 0.81321)*., 

represent unit vectors parallel to the direction of 
the strain and the normal to the habit plane 
respectively, and 

(6.5) m = 0.19144, 

is the magnitude (I-§5) of this strain. 

In their analysis of this strain, Greninger and 
Troiano assumed it to be a simple shear on the 
habit plane, and determined the shear direction 
as the intersection of the habit plane with a plane 
defined by direction of a specimen edge before and 
after transformation. If the theory is correct all 
such planes should intersect along the predicted 
direction (6.3). The observed scatter of Greninger 
and Troiano’s shear direction is just that expected, 
provided all specimen edges studied were within 
about 60° of the normal to the habit plane. The 
direction of their trial shear “A” is about 1° from 
the projection of the predicted direction on the 
habit plane. The measured angles of shear range 
from 8° to 14° with a mean of about 11° which 
compares favourably with the prediction that the 
direction normal to the habit plane rotates by 
11.0°. Thus, the theory is not inconsistent with 
these observations. 

7. Machlin and Cohen’s (259) -Transformation 

Machlin and Cohen’s [14] analysis of the trans- 
formation in Fe-30% Ni is the most recent 
measurement of the homogeneous strain accom- 
panying the transformation. The predictions of the 
theory for this case are summarized below using 
7 = 1 and 0 2 = 0.647. 

The invariant line strain for a and w positive 
now becomes 

/l, — 0.093129, 0.308397 \ 

(7.1) (lSl) = 0, 1.114712,0.016826 ), 

\0, 0.115080, 0.935472/ 

and cos a = 0.301977. The angle between [ll0] p 
and [111]b is now 3.1°, while the angle between 
(111) F and (101) B is 24'. This orientation relation- 
ship is about |° nearer the Nishiyama relationship 
than the previous one (§6). No precision deter- 
mination of the orientation relationship has been 
made but the predicted orientation does lie within 



the scatter of the results obtained by the pole 
figure method [6; 15]. 

The homogeneous strain is given by an expression 
of the form (6.2) where 

(7.2) - d = [0.23008, 0.63891, - 0.73407] P , 

(7.3) p' = (0.20226, 0.56166, 0.80227) F , 

(7.4) m = 0.22254, 

while the inhomogeneous shear is 0.4031 of the 
twinning shear. Machlin and Cohen estimate 
that the homogeneous strain has a shear compo- 
nent of 0.20 and a component normal to the habit 
plane of 0.05; the corresponding values predicted 
by the theory are 0.219 and 0.041. 

In comparing the predicted direction (7.2) with 
Machlin and Cohen’s observations, it must be 
remembered that they were only able to measure 
the projection of the direction on the (lll) F -plane. 
In Figure 8, their measured direction is compared 




Figure 8. Stereographic projection showing the predicted 
variation (for 7 = 1) of the direction of the homogeneous 
strain. The projection of this direction on the (lll)p-plane 
is compared with Machlin and Cohen's results. 

with the projection of the predicted direction on 
this plane. The direction given by Machlin and 
Cohen is a mean of all their results. Since the 
scatter of their individual measurements is con- 
servatively estimated, from the data in their 
Table I*, to be ±10° the agreement between 
theory and experiment is reasonable. The predicted 
direction is about 5° out of the (lll) F -plane which 
is consistent with the errors stated in Machlin 
and Cohen’s paper [16]. 

The scatter of the habit planes determined by 
Machlin and Cohen (Fig. 7) corresponds to the 
range of 0 2 from 0.64 to 0.66. The use of the value 

*The projection of d on the (lll) P -plane is [t X z/tyz, 1, 0] 
referred to axes in the directions [112] p , [110] P , [lll] p . 
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0 2 = 0.647 in the above comparisons is not of 
particular significance because the corresponding 
range of predicted directions lies almost entirely 
within the scatter of the experimental results 
(Fig. 8). If the scatter of the habit planes is real, 
it will be necessary to make measurements on 
individual martensite plates and then, perhaps, it 
may be possible to distinguish between the present 
hypothesis and those for which the complementary 
strain is either on the twinning plane or in the 
twinning direction. 

8. The (225) -Transformations 

In iron-carbon alloys containing less than 1.4% 
carbon the habit planes all contain a close-packed 
direction (110) F . Above 0.9%-carbon the habit 
planes are in the vicinity of (225) F (see Fig. 5), 
but as the carbon content is lowered below 0.9% 
the martensite plates degenerate first into laths 
and then into needles. Although the habit planes 
have not been determined in these cases the long 
axes of the crystals are in (110) F directions. 
Accordingly in applying the theory to these 
transformations d 2 is taken to have its maximum 
value. Thus, a = 90° and the invariant line is 
parallel to the twinning direction [110] F or [llI] B . 
Equation (4.10) for (lSl) now simplifies somewhat 
and it can readily be shown that the (a+, «+)- 
and (a+, w — )-solutions represent the same strain; 
similarly, for the other pair. Remembering the 
relations between solutions described in section 4, 
it follows that the variants are twin related in 
pairs. Hence, this group of transformations is a 
degenerate case of the general theory. 

Since [ll0] F and [111] B are now always parallel, 
the orientation relationships can be described 
conveniently by the angle \p between the planes 
(lll) F and (101 )b- The angle \p increases as the 
pole (101) B moves from (111) F towards (001) F 
and is given by 

( 8 . 1 ) $ = a + £ - n, 

where co is the angle (equation (4.11)) between 
(110) F and (112) b , i; the angle between (101) B 
and (112) b , and y the angle between (111) F and 
(110) F . Thus, 

sin* = ( 2 - T 2 )*[( 2 + 7 2 ) i -V 27 ]/ 2 ( 2 +y)n 

(8.2) sin f = 7/V2(l + 7 2 )*, 

sin i ? = 1/V3, J 

and the calculated variation of ip with y is shown 
in Figure 9. The predicted orientation relationship 
is exactly that of Kurdjumow-Sachs for 7=1 and 



never differs from it by more than about 15'. 
This is in agreement with the pole figure determina- 
tions of Kurdjumow-Sachs [3] and of Wassermann 

[17]. 




Figure 9. Comparison of the predictions of the present 
theory with that of Frank. $ is the angle between the ha bt 
plane and (111) F , and <p the angle between the planes (111) F 
and (101) B . 

The homogeneous strain is given by an expression 
of the form (6.2) where 

(8.3) -d = [(2-7 2 ) i ,(2-7 2 ) i , -27] F /V2(2+7 2 )* 

s [112] P , for 7 = 1, 

(8.4) p' = ((2-7 2 )l,(2-7 2 ) i ,V2(2+7 2 )l) F /2V2, 

- (2, 2, 4.90) F , for 7 “ 1, 

(8.5) m = 2[(2 + 7 2 ) 4 - V27]/(2 + 7 2 ). 

The inhomogeneous shear is always one-half of 
the twinning shear. It follows that twin orienta- 
tions with same habit plane are produced equiva- 
lently and are therefore variants as stated above. 

The above predictions for this special case are 
exactly the same as those obtained by Bowles [5] 
but differ somewhat from those of Frank [18], 
Frank suggested that the habit plane is determined 
by the condition that the close-packed planes 
(111) F , (101) B meet edge to edge in the habit 
plane. This is automatically true in the present 
theory provided that the (lll) F -plane intersects 
the habit plane (and the twinning plane) along 
the invariant line [110] F . This condition is satisfied 
in the present case (a = 90°) but will not always 
hold, so that Frank’s theory cannot be extended 
to more general cases. 

The predictions of the two theories are compared 
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in Figure 9. The angle <p (equation (8.1)) between 
the planes (111)f and (101) B is defined in the same 
way as in Frank’s paper, while his angle <t> between 

(111) F and the normal to the habit plane is here 
denoted by $ and, on the present theory, is 
given by 

(8.6) cos $ = [V2(2 - t 2 )1 + (2 + y 2 ) i ]/2 v / 3. 

The curves labelled “Frank” have been computed 
from the equations in §-7 of Frank’s paper using 
the uncorrected parameters given therein. The 
habit planes so predicted are clearly incompatible 
with observation. To improve the agreement 
Frank has suggested that before applying the theory 
the experimental lattice parameters should be 
adjusted so that the atoms in the [110] F , [llI] B 
rows match. One way of doing this is to change all 
the parameters of one structure by the same 
factor; this is equivalent to the introduction 
of the dilatation 8 in the present theory. When 
this is done the two theories give the same result. 
Frank employed a different method of correction 
but apparently applied the correction in the 
wrong sense for his procedure leads to even worse 
agreement with observation. 

In the development of his theory Frank has 
taken as clues the parallelism of the close-packed 
directions [110] r , [llI] B , and the close parallelism 
of the close-packed planes (111) F , (101) B . He 
envisages a simple type of interface in which these 
planes meet edge to edge along these directions. 
In the development of the present theory the 
authors have argued (Part I) that twinning is of 
fundamental significance. The twinning plane 

(112) B is generated from (I10) F and since, in the 
present theory, these planes always intersect 
along the invariant line they necessarily meet edge 
to edge in the interface in all cases. The agreement 
between the two theories for (225) F -transformations 
stems from the coincidence that in these trans- 
formations the invariant line is [ll0] F which also 
lies in the plane (111) P . 

9. Discussion 

The general conclusion to be drawn from the 
preceding analysis of specific transformations is 
that the theory is in substantial agreement with 
the available experimental data, although this 
data is inadequate for a conclusive test. To 
distinguish between the present hypothesis and 
the alternative ones where the complementary 
shear is either on the twinning plane or in the 
twinning direction, measurements of properties 



such as the habit plane, orientation relationship, 
direction and magnitude of the homogeneous 
strain, will have to be made on individual marten- 
site plates. 

The present hypothesis provides a description 
of the overall atomic displacements produced by a 
transformation, which is certainly very nearly 
correct. If it is correct, the problem of finding 
suitable mechanisms for these transformations 
is reduced to discovering those which are compat- 
ible with the hypothesis; any such mechanism 
would automatically be consistent with all the 
geometrical properties of the transformations. 
It seems almost certain that a satisfactory mech- 
anism will involve the motion of dislocations and 
the theory suggests that these dislocations will be 
related in some way to twinning in the final 
structure. 

The theory contains one adjustable parameter, 
8, which has been estimated from the observed 
habit plane. It may be asked what physical or 
geometrical requirement determines its value. 
If the requirement were purely geometrical, it 
might be hoped that its nature would emerge 
from a study of the transformations as a whole. 
Although this has not proved to be the case, an 
interesting feature of the iron -carbon group of 
transformations has been noticed and may be of 
some significance. 

The habit planes observed in iron-carbon alloys 
containing 0-1.78% carbon are plotted in Figure 
10, together with the predicted variation of the 
habit plane, for 0 = 54° 44' (cos /3 = 1/V3), and 
0 = 56° 20'.* Along these latitude circles the 




Figure 10. The same as figure 5, together with the lines 
along which 0 takes the values 54° 44' and 56° 20'. 



•For a = 90°, tan (54° 44' — <f) — cos 0, and then these 
values of 0 correspond to habit planes 30° and 29° from 
(001) F respectively. 
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angle ^ between the close-packed planes is zero 
and approximately 7' respectively. Although the 
data is inadequate, it is in agreement with a 
requirement that as the carbon content increases 
\f/ increases to about 7' and thereafter remains 
approximately constant. As y increases the trans- 
formations can be of the (225) P -type until \f/ 
reaches 7' ; thereafter, the habit plane would move 
rapidly out into the centre of the triangle along 
the line /S = 56° 20', for the lines y = const, 
cross the lines /3 = const, at a small angle. Such a 
rapid change of habit plane with y could account 
for the reported discontinuous variation of the 
habit plane at 1.4% carbon [19]. 

Another feature of interest is that the latitude 
circle corresponding to cos /3 = 1 / \/3 and \p = 0, 
passes through (011) F - If this latitude circle has 
any particular significance one might expect 
that a transformation with y about y/2 would 
have a habit plane in the neighbourhood of (011) F . 
This is the case for transformations such as In-Tl 
[20; 21] or Cu-Mn [22]. These transformations are 
special cases of the present class. 
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